TRAVELING WAVES IN A THREE SPECIES 
COMPETITION-COOPERATION SYSTEM 
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Abstract. This paper studies the travehng wave solutions to a three species 
competition cooperation system. The existence of the traveling weaves is inves- 
tigated via monotone iteration method. The upper and lower solutions come 
from either the waves of KPP equation or those of certain Lotka Volterra 
system. We also derive the asymptotics and uniqueness of the wave solutions. 
The results are then applied to a Lotka Volterra system with spatially averaged 
and temporally delayed competition. 



1. INTRODUCTION 

We study the traveling wave solutions of the following three species competition 
cooperation system 

{Ut = Uxx + u(l - u - aiw), 
Ut=Vxx + rv{l~a2U-v), (x,t)eRxR+. 
Wt = Wxx + -{v - W), 

where u{x,t), v(x,t) and w(x,t) stand for the population densities of the three 
different species, > is interaction constant, i — 1,2 and r > (— ^ < 0, 
respectively) is the relative intrinsic growth rate of the species v (w, respectively). 
Aside from the intra-specific competitions, system (|l.ip describes the relation that 
the species w competes with u and u competes with v, while v cooperates with w. 

The purpose of our study is of two folds: there are less results (see[T7l |9]) on 
the traveling wave solutions to the three species systems even for one with simple 
form as (jl.ip ; on the other hand we would like to extend the results of [T] from the 
tempo-spatial delayed KPP (Kolmogorov, Petrovsky and Piscounov) equation to 
the following Lotka Volterra competition system with spatial temporal delay 



(1-2) 



Ut = Uxx + u{l ~ u - aig * *v), 
Vt = Vxx + rv{l - a2U - v), 



where the function n * *v — r^°° , e -e * ^ ' v(s)dsdv repre- 

y J-oo J-oo ^2TT(t-s) ^ \ I y f 

sents the tempo-spatial delay of the response of species v to u, here g{x,t) = 
—i — eTt-e~i, and satisfies 



Key words and phrases. Traveling wave, spatio-temporal delay, Lotka Volterra, Competition 
Cooperation, Existence. 

1 



WAVES IN THREE SPECIES SYSTEM 



2 



(1.3) g * *1 ~ 1. 

Noting that ii g * *v is replaced by v in (|1.2|) . we recover the classical Lotka 
Volterra competition system, and fruitful results have been devoted in the study of 
traveling waves [3l[5l|8l[2l[l3l[l4l[TTl[12l[15l[T9l[2Q] arising from it. It is interesting 
to see the long term effect of introducing spatio-temporal delay to the competition. 
The temporal delay accounts for the time once consumed resource the dominated 
species needs to wait for its re-growth, and the spatial averaging accounts for the 
fact that individuals are moving around and have therefore not been at the same 
point in space at different times in their history (see [I]). 

On setting 

1 /■ + °° 1 (x-i,)^ _(t-s) 

(1.4) 'w{x,t) = - / / _ e4(t~s)e . v{s)dsdy, 



T 



-oo y/2TT{t — s) 



we easily verify that w{x,t) satisfies the following equation 

(1.5) wt ^ Wxx + -{v - w). 

T 

System (|1.2p is now recasted into (II. ip . The existence of the traveling wave 
solution of (II. ip is equivalent to that of (II. 2p . 

Throughout the paper we make the following assumptions: 

(1.6) Condition HI : < ai < 1 < 02 
and either 



(1.7) Condition H2a : r{a2 - 1) < (1 - ai) 



(1.8) Condition H2b : r(a2 - 1) > (1 - ai) > r{aia2 - 1) 



(1.9) Condition H3 : 02 > is suitably large . 

Under the condition ()1.6p . system (jl.ip admits three constant equilibria: (0, 0, 0), 
(0, 1, 1) and (1, 0, 0), and the first two are unstable while the third is stable. We are 
interested in finding the traveling wave solutions of (jl.ip connecting (0, 1, 1) with 
(1,0,0). Via transformation (|1.4I) . the existence of the traveling was solutions of 
(jl.ip connecting (0,1,1) to (1,0,0) is equivalent to that of (|1.2p connecting (0,1) 
with (1,0). 

A traveling wave solution for (jl.ip has the form (w(x, t), v{x, t), w{x, t)) = {u{x-\- 
ct), v(x + ct),w(x + ct)) = (w(0, w(0: — X + ct^ and satisfies the system 

— cu^ -\- u{\ — u — aiw) = 0; 

— cv^ + rv{l — ^ v) — 0, 
(1-10) { 1/ ^ n 

' w^^ ~ cw^ H — (V ~ w) ~ V, 

{u,v,w){—(X)) — (0,1,1), {w,v,w){+(X)) — (1,0,0). 
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For the convenience of later study, we change (|1.10|) into monotone ( p3] ) . Let 
u^u, v = l — V and w — I — w, and drop the bars on u, v and w, we arive at 

— cu^ + u{l — fli — u + aiw) — 0, 
v^^ — cv^ + r{l — v){a2U — u) = 0, 

(1-11) { 1/ N ^ 

w^^ — cw^ H — [V — w) =0, 

{u,v,w){-oo) = (0,0,0), {u,v,w){+oo) = (1,1,1). 

Conditions (|1.7I) and (|1.8p stipulate that wave solution is either below or above 
the plane u = v, and this is reflected in the construction of the upper and lower 
solutions in section [31 

Definition 1. A C'^{R) x C'^{R) function {u{^),v{£,),w {£,)), C e K is an upper 
solution of p. lip if it satisfies the inequalities 



(1.12) <^ 



ujj — cu^ + u{l — ai — u + aiw) < 0, 
v^^ — cv^ + r{l — v){a2U — v) < 0, 

w^^ — cw^ H — (v — w) < 0, 





T 

and the boundary conditions 

(1.13) 1^ V j (-00) > ^ I , ( V I (+(x)) > 

A lower solution of (jl.lip is defined in a similar way by reversing the inequalities 
in ([1121) and ([1131) . 

Remark 2. The minimum of two upper solutions is still an upper solution, i.e. 
suppose Wi(^) = (wi, W2, W3)(^) and W{^) = (wi, W2, W3)(^), ^ € M are two upper 

solutions, then V{^) = (min(wi(^), wi (^)), min(?«2(0, u)2(0): niin(w'3(C): ^^3(0)) is 
also an upper solution. For the lower solutions, we have similar observation except 
that the minimum is replaced by maximum. 

Since (jl.lip is a monotone system, the monotone iteration method ([2Q]) is ready 
to apply once the orderness of the upper and lower solutions is estabished. The key 
to the monotone iteration is to identify a pair of ordered upper and lower solutions 
to (|l.lip ([20l|3]). There are two methods to set up the upper and lower solutions. 
The first one was used in [3l [3 [20] , which consists of a pair non-smooth upper and 
lower solutions, and the similar idea was later successfully generalized to handle 
local and nonlocal equations, the second one is based on a pair of smooth upper 
and lower solutions from known equations, and the method was applied in |15j for 
a general form of two species Lotka Volterra competition system and in [6] for a 
model system arising from game theory. We will use the ideas of the second method 
to set up the upper and lower solutions for (|l.lip . See section 2 for details. 

The other interesting aspects of the traveling front solutions are the minimal 
wave speed, the uniqueness, the asymptotics and the stability. The minimal wave 
speed is also referred to as the critical wave speed, below which the there will be no 
monotonic traveling waves. Also the traveling waves with the critical speed behaves 
differently at —00, see [15) . We will use a generalized version of sliding domain 
method (see [2] ) to show the uniqueness of the front solution corresponding to 
each speed. 



WAVES IN THREE SPECIES SYSTEM 



4 



This paper is organized as follows: In section [2] we gather the necessary in- 
formation about the KPP and Lotka Volterra waves, in particular we derive the 
asymptotics of the Lotka Volterra waves at —oo which is the key in setting up the 
upper and lower solutions for In section [3] we show the existence of the wave 

solutions for (ll.lip and further derive their properties such as strict monotonicity, 
the uniqueness and the asymptotics. 

2. PROPERTIES OF WAVES FOR KPP AND A TWO SPECIES LOTKA VOLTERRA 

COMPETITION EQUATIONS 

In this section, we introduce porperties of the wave solutions to KPP equation 
and to a two species Lotka Volterra competition system, which will be a key ingre- 
dient in the construction of the upper and lower solutions for system (|l.lll) . For 
the rest of the paper the inequality between two vectors is component-wise. 

The construction of the smooth upper and lower solution pairs for system 
is based on the known results on the KPP equations and the recent results on a 
rescaled Lotka Volterra system. It seems that the asymptotics of the Lotka Volterra 
waves derived in this section is new. 

Consider the following form of the KPP equation: 

r uj" ~cuj' + f{uj) = 0, 

(2.1) <^ 

[ w(— oo) = 0, a;(-foo) = 6, 

where / e C2([0, h]) and / > on the open interval (0,6) with /(O) = /(&) = 0, 
/'(O) = a > and f'{b) = -bi < 0. We first recaU the following result (JTE\): 

Lemma 3. Corresponding to every c > 2-\/^, system 12. 1\) has a unique (up to 
a translation of the origin) monotonically increasing traveling wave solution uj{^) 
for ^ € M. The traveling wave solution w(^), ^ € R has the following asymptotic 
behaviors: 

For the wave solution with non-critical speed c > 2\/d., we have 

(2.2) w(^) = a^e" + o(e' ^— €) as ^ ^ -oo, 

c--\/c-^+4fji ^ c-\/c^-|-4fai - 

(2.3) iu{0 = b - b^e ^ « + o(e ^ «) as ^ +oo, 

where a^j and b^i are positive constants. 

For the wave with critical speed c — 1\f^, we have 

(2.4) w(0 = (flc + dcOe^^ + o(Ce^«) asi^ -oo, 

(2.5) a;(0 = b - fo^e^^-^^'? + o{e^^- '^^^^^) as i ^ -^oo, 

where the constant dc is negative, be is positive, and Oc G M. 

We also need the existence and asymptotics (at — oo) of the solutions of the 
following rescaled version of Lotka Volterra system: 
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(2.6) 



u^^ — cuj + u(l — ai — u + aiw) — 0, 
v^^ — cv^ + r(l — v){a2U — v) = 0, 
^ (m,i;)(-oo) = (0,0), (u,«)(+c^) = (1,1). 



The asymptotics of the wave solutions wih be derived by comparing the as- 
ymptotic decay rates of the upper and lower sohitions of (|2.6p . Recalhng that 
(w,w)(0 € C2(M) X C2(R), ^ e R is an upper solution for dH]) if it satisfies 



u^^ — cu(^ + u{l — fli — w + aiw) < 0, 

v^^ — cv^ + r(l — v){a2U — v) < 0, 

(w,i>)(-(X3) > (0,0), (u,t;)(+oo) > (1,1). 

by reversing the above inequahties. 



(2.7) 



and we can define the lower sohition for 
We have the fohowing conclusions: 

Lemma 4. Let the parameters satisfy either HI, H2a or HI, H2h, then for each 
c > 1\J\ — a\ , system \2. 6\) has a unique (up to a shift of the origin) strictly 
monotonic solution, and for < c < 2-^/1 — ai, 112. 6]) does not have monotonic 
solution. At — oo the solution has the following asymptotical properties: 
For c > 2^1 — fli , the solution {u{^),v{^)) satisfies, as ^ ^ ~oo; 

"(C) \ Ml 



(2.8) 



+ o(e" 



Vc2-4(l-ai) , 



while for c* = 2\/l — ai, the solution (u(^), w(^))"^ satisfies 
"(0 \ / Anc + Ai2c^ 



(2.9) 



o(Ce 



V{0 J \ A21c + A22c^ 

as ^ — > — oo, where Ai, A2 > 0, Anc, A21C G R and Ai2c, ^22c < 0. 

Proof. The existence of the waves under conditions Hl-H2a, or HI, H2b is contained 
in [TB] , and also refer to [T3] for the existence and the asymptotics of the waves in 
a more general form of a two species competition system under conditions Hl-H2a. 

From now on we will concentrate on system (|2.6p under conditions HI and H2b. 
The authors in [16] proved the existence of the monotone solutions for c > 2y/l — oi 
by showing that (|2.6|) is linearly determinated. However such method does not bring 



us the crucial information on the aymptotics of the wave solutions that is needed in 
section [31 We show that the traveling wave as derived in [TB] is actually squeezed 
by the lower and upper solutions of (12.61) constructed below. Noting the upper and 
lower solutions differ from that in j8]. 

We first set up the lower solution for (12.61) . 

For a fixed c > 2^yi — ai, let m(^), ^ € R be a corresponding solution of the 
following KPP equation 
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(2.10) 



u" - cu' + (1 - ai)u(l - u) = 0, 
u{~Qo) — 0, u{+oo) = 1. 



It is straightforward to verify the fohowing claim. 



Claim A. Under conditions HI and H2b, for each fixed c > 2^/1 — oi, the smooth 
function (w, = {UjlLjiS,)^ ^ defines a lower solution for Ii2.6\) . 

We next set up an super sohition for (|2.6p . Choosing a smah number I such that 



(2.11) 



0<l< 



1 - ai - r{aia2 ~ 1) 



1 + r — ai 

Let c > 2V1 — fli be fixed and m(^) be a solution of the following modified KPP 
equation; 

I 



(2.12) 



w" - cu' + (1 - ai)w(l - 
it(— oo) = 0, u(+oo) = 



1 — fli 

1 — ai 



-u) = 0, 



I 



Condition H2b (see II. 8p implies that 



1 — ai — r(aia2 — 1) 
1+r-ai 



< 1 — fli we therefore have 



fi>l. 
Setting 



(2.13) 



(^(0,^«(0), ifu(0<T^; 

a\ 1 — t 

(u(C),i), ifY^<zi(0<i; 



I (1,1), if^i(e)>i. 

Claim B . Assume conditions in Claim A and let I satisfy k2.11\) . then k2.13\) defines 
an upper solution for \2.6\) . 

Proof of claim B. It is easy to see that {u,v)(—oo) = (0,0), {u,v){+(X)) = (1,1), 
and (u,u)(^) = (1, 1), ^ € M solves the first two equations of (|2.6I) . We next verify 

that {u,v){^) — (u(^), w(C)) is also an upper solution for (|2.6p . 

ai 

For the u component we have 

m" — cu' + u{l — fli — u + aiv) 

— ii" — cu! + m(1 — ai — u + (1 — l)u) 
= u" ~ cii' + {I - ai)u{l - Y^u) 

— (1 — ai)u{l — j-t^m) + ■"(I — fli — u + aiw) 
= m[1 - ai - M + (1 - 0" - (1 - ai)(l - T^^u)] = 0, 
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and for the v component we have 

f)^^ — cv^ + r(l — v){a2U — v) 

= hzlu'' ^ c^^u' + r(l ^ ^^u)(a2U - ^u) 



(2.14) 

+ (1 - ai)u(l - j^u) - (1 - aiHl - j^u)] 
- liriaiKl - lir")(«2 - - (1 - ai)^(l - j^u)] 

= IIT^IKt^ - 1) - (1 - ai) + HI ~ r{a2 - i^))]. 
To ensure the non-positiveness of the last expression m (j2.14l) . we require: 

(2.15) ,(^_i)_(i_„,)<0, 
and either 

(2.16) l-r{a2-^—^)<0, 

ai 

or 

l- r(a2 ~ ^) > 0, and 
r{fzf - 1) - (1 - ai) + iiiOil ~ r{a2 ~ ^)) < for ah ^ G 1 



It is a tedious but straightforward verification that for I satisfy (|2.1ip . we have 
(|2.15p and (|2.16p . or (I2.15P and (|2.17p . Hence in either case the last expression in 
(|2.14p is less than or equal to 0. 

Analogous to Remark[2]we see that (|2.13p is an upper solution. This finishes the 
proof of the Claim B. 

Noting the relation 

(2.18) uiO = ^-f^i^iO > uiO e e K, 

there follows the orderness of the upper and lower solution 

(2.19) (uio.viO) < eeM. 

By the monotone iteration scheme ([20]), for each fixed c > 2-^/1 — ai the traveling 
wave solution (u,w)(^) satisfies 

(2.20) {u,vm < {n,vm < ? e R. 



Let c > 2\/l — ai be fixed and let {u{^),v{^)), ^ G M be a corresponding solution. 
The upper and lower solutions have the exactly asymptotic decay rate at — ooby 
Lemma [3l the estimates (|2.8p and (|2.9p then readily follow. □ 



3. EXISTENCE OF THE WAVES 



3.1. Ordered upper and lower solutions under conditions HI and H2a. 

To construct the upper-solution for the system (|l.lip in this case, we begin with 
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the following form of KPP system: 

( u" - cu' + (1 - ai){t(l - u) = 0, 
(3.1) <^ 

[ oo) = 0, u(+oo) = 1, 

where relating to dH]), /({)) = (1 - ai)w(l - u) > for u G (0, 1). /(O) = /(I) = 0, 
/'(O) = (1 — fli) > and /'(I) = —(1 — ai) < 0. According to Lemma|3l for each 
fixed c > 2^/1 — ai, system has a unique (up to a translation of the origin) 
traveling wave solution u{S,) satisfying the given boundary conditions. Define 



(3.2) 



V m j 

we have the following result, 



Lemma 5. Assume the conditions HI and H2a, for each fixed c > 2^1 — ai, i3. 2(l 
is a smooth upper solution for system 

Proof. On the boundary, one has (u, w, i(})(— oo) = (0,0,0), {u, v, wY' {+oo) = 
(1,1,1). 

As for the u component, wc have 

u" — cu' + u{\ — ai — u + aiw) 
= u" — cu' + u(l — ai — u + aiu) 
= u" — cu' + (1 — ai)u(l — u) 
- 0, 

and for the v component 

v" — cv' + r(l — v){a2U — v) 
— u" — cu' + r(l — u){a2U — u) 
+ (1 — ai)u{\ — u) — (1 — ai)u(l — u) 
= [r(a2 - 1) - (1 - ai)]?2(l - u) < 



due to the condition H2a. 
As for the w component. 



w — cw 



v' + ^{v-w) 
u" — cu! + (1 — ai)u{\ — u) — (1 — ai)f2(l — u) 



= -(1 - ai)f2(l -m) < 0. 
Thus the conclusion follows. 



□ 
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We next construct the lower solution pair for system For any small but 

fixed number I with 



(3.3) < / < 



we choose a number I such that 



ra2 



1 — oi + r ' 



(3.4) 



< ; < 



I 



(l-ai)r + l' 
We begin with yet another KPP system: 



(3.5) 



ii" - cu' + (1 - ai)M(l - u) = 



u{—oo) — 0, u{+oo) 
Corresponding to the notions in Lemma [3l 



1 — ai 

1 — fli 
l + l 



< 1. 



f{u) = (1 - ai)u(l - ^-^u) > 
1 — ai 

for u e (0, ^-^). /(O) = /(^^) = 0, /'(O) = 1 - ai > 0, and /'(\^) < 0. 
For each c > 2^1 ~ ai let w(^), ^ e M be a solution of (|3.5p and let 

/ m \ I u{i) \ 



(3.6) 



we have 



V HO J 



Lemma 6. For each c > 2^1 — ai , \3. 6]) is a smooth lower solution of system 



Proof. On the boundary, one has 



[ u(-oo) ^ 
v_{—oo) 



( \ 





/ 2i(+oo) \ 

u(+oo) 
y w(+oo) ) 



1 + / 
1 — ai 

TTT 

1 — fli 



< 



/ 1 \ 

1 

V 1 / 
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Furthermore, for the u component, 

u" — cu' + u(l — ai — M + aiw) 

= u" — cy! + u(l — fli — M + ailu) 

+ (1 - ai)M(l - u) - (1 - ai)u(l - u) 

1 — ai 1 — ai 

= — fll — u + fliZu — (1 — fll) + (1 + 

> 0, 

and for the v component, we have 

— cu' + r(l — v) {a2U — v) 

T 

= — cu' + y (1 — lu){a2U — lu) 

+ (1 - ai)u(l - — — m) — (1 — ai)u(l - ^-^u)] 
1 — Oi 1 — ai 

I 1 — ai 

= lu{ja2 - r - (1 - ai) + [(1 + I) ^ r{a2 - l)]u} > 

because of conditfon p.Sp . 

As for the w component we have 

w" — cw' -i — {v — w) 

T 

= Iv/' — clu' H — {lu — lu) 

T 

= lu[-Al^l)-il-ai){l-^^u)] 
= lu[\{l-l)-{l-ai) + {l + l)u]>0, 

Tl 

due to the choice of I and I. 

The conclusion of the lemma follows. □ 

Remark!. 1. For any fixed c > 2y/l — ai and {u, v, 'w)'^{^), {u,v,w)'^{£,), ^ £ M the 
respectively upper and lower solutions defined in p.2p and p.4p . we have relation 

{u,v,w)'^{0 > {u,v,w)'^{0 for C e R. In fact u{0 = W^(^) < "(^)' ^ ^ 
where u(^) and are the solutions of (|3.5p and p.ip respectively, it then follows 

m = ^1^(0 < s(o = m and = HO < m = ^io- 
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2. The construction of lower solution in Lemma 5 also applies to the case Hl- 

TCi'2 

H2b, where we simply require the condition < Z < to be replaced by 

\ — ai+ r 

rar> 

0<l< min{ , 1}. 

1 - ai + r 

3.2. Ordered upper and lower solutions under conditions HI, H2b and 

H3. The following estimates of the solutions of system (|2.6p is needed in the con- 
struction of the upper and lower solutions. 

Lemma 8. Assume the conditions of Lemma ^ Let (u(^),w(^)), ^ G M, be a 
solution of \2.6\) for a fixed c > 2\/l — ai, then there exists a constant 02 > such 
that for 02 > a2u{^) > v{^) for ^ £ R. 

Proof. Noting that u{£^) and w(^) have the exactly same (up to the first order) 
exponential decay rate at —00 , and at +00 we have a2u{^) > v{^). The rest of the 
proof follows easily. □ 

We now set the upper and lower solution pairs for system 



Lemma 9. Let the parameters satisfy HI U.6\) and H2b U.8]} . then 113. 6]) consists 
of a lower solution for iLll]) . 

Proof. Similar to the proof of Lemma |6] so we skip it. □ 
We next set up the upper solution for 

For each fixed c > 2^/1 — ai let {u,v){£,) be the correspondingly the unique 
solution of p.41|) . and we define the function 

(3.7) iu,v,wm^iu,v,vm, 
We then have 

Lemma 10. Let the parameters satisfy HI U.6]) . H2b il.8]) and H3 !il.9\) . and 02 
be given in Lemma\^ then for a2 > ftj ^3. 7[ ) defines an upper solution of 

Proof. The proof follows easily from the fact that {u,v){£^), ^ e M solves (|2.6p and 
Lemma [S] □ 

We next show that such constructed upper and lower solutions are ordered. The 
following generalized version of sliding domain method ([!]) is needed. 

Proposition 11. Let two vector functions U{£,) = (ui{£^),U2{£,), ...jUniO) '^'^'^ 
IZ(C) — (l*i(?):l*2(0i (0) satisfy the following inequalities: 

(3.8) DU" - cU' + F{U) < < DU!' - c]f + F([/) for ^ e [-N, N] 
and 

(3.9) m-N)<UiO forU{~N,N], 

(3.10) mo < U{N) for C e [-N, N), 

where D is a diagonal matrix with positive entries Di, i — 1, 2...n, F{U) = {Fi(U), Fn 
is with respect to its components and > for i 7^ j , i,j — 1, 2...n, then 
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(3.11) mO<U{0, (e[~N,N]. 

Proof. We_ adapt the proof of 2 . Shift U{^) to the left, for < ^ < 2N, consider 
U'^iO ■= U{£, + A*) on the interval {—N — ji^N — ji). On both ends of the interval, 
by (pr^ and (PIIT)) . we have 



(3.12) 



Starting from fi — 2N, decreasing /i, for every /i in < < 2N, the inequality 
p.l2p is true on the end points of the respective interval. On decreasing /i, suppose 
that there is a first fi with < /i < 2N such that 

and there is one component, for example the i — th, such that the equality holds on 
a point ^1 inside the interval. Let W{S,) = (wi(^), W2(0j •■•7 '"^n(O) = U^{^) — IZ(C)j 
then Wi{£,), i — 1,2, n satisfies 



ii;.(a)-0, w,{^)>Q iorie[-N-^i,N-^jL], 
the Maximum principle further implies that = for ^ G [—N — ^l, iV — /i], but this 
is in contradiction with p.l2p on the boundary points ^ = —N — /i and ^ = iV — /i. 
So we can decrease /j, all the way to zero. This proves the conclusion. □ 

We next shift the upper solution obtained in Lemma [TU] far to the left to achieve 
the orderness between the upper and lower solutions. 



Lemma 12. Let c > 2^/^~^^ll he fixed and {u,v,w){S,), {u,2Li2il){C) the corre- 
sponding upper and lower solutions derived in Lemma \T(]\ and Lemma\^ then there 
exists a r/o > such that for all rj > rjQ we have 

iu,v,w){^ + f])>iu,v,w){^), CeK. 

Proof On the boundary, we have (m,w,w)(^) (1, 1, 1) and (u, u, w;)(C) {^Y^^^^Tft^^^TTt) < 
(1,1,1) as ^ -> +00 , hence there exists a suSiciently large A^i > such that for 
any t] > 0, 



(3.13) iu,v,w){^ + v)> iu,21,wm, ee[iVi,+oo). 

while at ^ = —00, we have the asymptotics of the upper and lower solutions: 
For c > 2^1 - ai. 



(3.14) 



and 



A2 

V A, j 



-V'c^-4(l-ai) , 



+ o(e" 



-Vc^-4(l-ai) 
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(3.15) v{0 

V iRiO J 

While for c = 2^/1 - ai 



(3.16) 



and 



(3.17) 



V m J 



( All + Sue \ 

^21+^21^ 
\ ^31 + ^31$ j 

( All + Sue \ 
A21 + B2ii 



« + o(e 5 «); 



gvTT^e + 0(6^"^^^) 



\ A31 + SaiC / 

Then it is easy to see that there exists a 770 > such that for any t] > rjo, we 
have a — iV2 < and the relation 

(3.18) iu,v,iv)i^ + Tj) > {u,v,wm, ^e{-(x,-N2] 

holds for c > 2vT~~a]". We may further adjust A^i and N2 such that p.l3p and 
p.lSp hold on the interval [N, +00) and (—00, N] respectively for some large > 0. 
While on the interval [—N, N], {u, v, w){(, + ri) and {u,v,w){^) satisfy the conditions 
of Proposition [m we therefore have 

(3.19) {u,v,w)i^ + v)>{u,v,wm, U[-N,N]. 

The conclusion of the lemma follows from (13.131) . p.lSp and (13. 19^ . □ 

In the sequal we still write the shifted upper solution as {u,v,'w){£,), ^ e R. 

3.3. Monotone waves and their asymptotics. With such constructed ordered 
upper and lower solution pairs, we now have 

Theorem 13. Assume either the conditions HI and H2a or HI, H2h and H3, then 
for every c > 2^/1 — ai, system has a unique (up to a translation of the 

origin) traveling wave solution. The traveling wave solution is strictly increasing 
on M and has the following asymptotic properties: 
1. Corresponding to the wave speed c > 2^/1 — a\, 

( <0 \ f A,\ 

c-v'c^-4(l-ai) 



(3.20) 



as e 



A, 



^ + o(e 5 ?) 



V w{0 J \A,J 



-00; 
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while Corresponding to the wave speed c = 2-y/l — ai, we have 



(3.21) 



as ^ 



( "(0 \ 



( (Aiic+Ai2cO \ 



For any speed c > 2^1 — ai, we /lawe 



(3.22) 



(3.23) 

M(0 



( 1 \ 



/ 1 \ 

1 

V 1 ) 



Sl 



e 2 



S2 



V / 
/ 1 \ 



i-(l-a2) + l \ 

rr(l - 02) + 1 

1 y 



Ve2+4(e.2-l) , 



+ S3 



V 1 / 



and the real numbers Si, i = 1,2,3 are not all zeros and be such that — M(^) + 
o(A/(^)) < for ^ > large enough. The speed c = 2^/1 — oi is the minimal wave 
in the sense that below which there is no monotone waves of 



Proof. Starting from the upper and lower solution pairs obtained in section \3A] and 
section 13.21 and using the monotone iteration scheme provided in j20[ |3] , we obtain 
the existence of the solution (u(^), -y(^), ?i'(^)) to (|1.11|) for every fixed c > 2^/1 — ai. 
The solution satisfies 



V MO J 



< 



V MO J 



< 



V MO J 



Lemma[3]and Lemma [4] imply that the upper- and the lower-solutions as derived 
in section [Ql and section [3T2] have the same asymptotic rates at —00. Then p.20p 
and (|3.2ip then follow from Lemmas [3] and ID 
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To derive the asymptotic decay rate of the travehng wave solutions at +00, we 
let c > 2^1 — fli and 

(3.24) C/(0 :=K0,«(0,«^(0) 



be the corresponding traveling wave solution of (jl.lll) generated from the mono- 
tone iteration. We differentiate (|l.lip with respect to and note that U'{(,) := 
{wi,W2,W3)'^ {^) satisfies 

(3.25) (wi)^^ - c{wi)^ + Aii{u,v,w)'Wi + Ai2{u, v, w)w2 + ^13(^^1 v, w)w3 = 0, 



(3.26) (W2)« - c{w2)^ + A2i{u,v,w)wi + A22{U,V,W)W2 + A23{U,V,W)W3 = 0, 



(3.27) (^3)?? - c(w3)^ + A3i{u,v, 10)11)1 + A32(m, U, W)W2 + A33{U,V,W)W3 = 
where 



All 


A12 


A 


A21 


A22 


A 


A31 


A32 


A 



13 

23 I = 

33 



f — fli — 2m + aiw 
a2r{l - v) 




aiu 
-r — a2ru + 2rv 



The limit system of ([3^ . ([3?26l) and ([3^ at ^ = +00 is 

(V'i)« - ciipi)^ - ipi + aiips = 0, 



(3.28) 



. (^3)« - C(V'3)« + 7^2 - 7^-3 = 0. 

It is easy to see that system p.28p admits exponential dichotomy (U)- Since 
the traveling wave solution (u(^), w(^), w(^)) converge monotonically to a con- 
stant limit as ^ — ±00, the derivative of the traveling wave solution satisfies 
(wi(±oo),u'2(±oo),it;3(±oo)) = (0,0,0) (^Q], p. 658 Lemma 3.2). Hence we are 
only interested in finding bounded solutions of p.28p at -l-oo. 

Introducing transformation ^/ = PY by 




1 


ai 

r(l-a2) + l 


aiT \ 


/ 2/1 (0 





rr(l — 02) + I 


Y 


2/2(0 





1 




\ 2/3(0 




(3.29) 

we can decouple (|3.28p into the following equivalent system: 

(2/1)?? - c{yi)^ -2/1=0, 
(3.30) I (2/2)?? -c(y2)e+r(l- 02)^2 =0, 

, (2/3)?? - ciys)^ - iyg = 0, 

and find its bounded solutions at -l-oo explicitely. In fact, for some nonzero con- 
stants di,d2,d3, we have 
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(3.31) yi(e) = die- 



2/2(0 = d2e~ 



2/3(0 = rfse" 



-VI 



Transforming back to wc have 



(3.32) 



MO 

MO \=P\ 2/2(0 I = 



MO 



2/i(0 \ / 2/i(0 + ViT^MO + 7^2/3(0 

(rr - Tra2 + 1)2/2(0 

2/3(0/ V 2/2(0 + 2/3(0 



Hence we have p.22p on intergrating p.32p . 

We next show the strict monotonicity of the travehng wave solutions, which 
wiU be a key ingredient in locating the eigenvalues of the linearized operator 
about the traveling wave in a separate study. By the monotone iteration pro- 
cess (see [20]), the traveling wave solution f/(0 is increasing for ^ e R, it then 
follows that iwi{^),W2iO,W3{OV = U'iO > and satisfies (jX^ . and 
(|3.27p . The monotonicity of system (ll.ll[) and the Maximum Principle imply that 
(wi, W2, W3)"'"(0 > for ^ e M. This concludes the strict monotonicity of the 
traveling wave solutions. 

On the uniqueness of the traveling wave solution for every c > 2-^/1 — ai, we 
only prove the conclusion for traveling wave solutions with asymptotic rates given in 
p. 201) and p.22p . since other case can be proved similarly. Let ?7i(0 = {ui, vi, wi)"^ 
and ?72(0 ~ ("2, V2, ^2)^ be two traveling wave solutions of system (|l.lip with the 
same speed c > 2y/l — ai. There exist positive constants Ai, Bi, i = 1,2,3,4 and 
a large number > such that for ^ < — A^, 

/ ^1 \ 

c+Vc^-4(l-ai) , 



(3.33) 



UiiO = 



(3.34) 
and for ^ > iV 
(3.35) 



c/2(0 



A2 

V ^3 J 

( A, \ 
A, 

V A, J 



« + o(e 5 ?) 



c+Vc^-4(l-ai) , 



+ o(e 



:;+\/c^-4(l-ai) , 



c/i(0 



1 - B2e^2« 



/ 1-Bie^^i« \ 

26 

/ 1-^1 



(3.36) 



t^2(0 = 



1 - B,e^2« 



V 1-^3 



/ o(e^i«) \ 
o(e^2«) 

V o(e^^3«) / 

o(e^2«) 

V o(e^3«) / 
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where /i^ is one of the elements in the set { "^^ ^^^^ , \/c^+4 ^ fL_y^EilZl}^ 
and Bi, are positive numbers, i — \^2, 3. The travehng wave solutions of system 
(11.111) are translation invariant, thus for any 9 > Q, C^i (f) + is also a 

traveling wave solution of (ll.lip . By p.33p and p.35p . the solution + 9) has 
the asymptotics 



(3.37) Ul{0 



for C < "N- 



V ^36 





2 




2 


■-^1 ' 




-4(1- 




-4(1- 


ai) 




2 




2 






-4(1- 




-4(1- 


ai) 



+ o(e 2 



(3.38) 



c/f(^) = 



/ 1 - Bie^i'e^i^ \ 
V 1 - B3e^3«e^3« y 



for C > A^- 

Choosing 6* > large enough such that 



(3.39) 



3+ij 



/ o(e^i«) \ 
o(e^2«) 

V o(eH3«) y 

i = l,2,3. 



(3.40) Bie^^^<B.^, 1,2,3 
then one has 

(3.41) C/f (0 > [72(C) 

for C S (— oo, — iV]U[Af + oo). We now consider system (|1.3p on [—N, +N]. We can 
verify that Uf{£,) = U{£,) and (C) = CZ.(C) satisfy all the conditions of Proposition 
[TTl hence we have f/f (C) > C^2(C) for C S [—N, N]. Further applying the Maximum 
Principle and noting that uf{±N) > C/2(±iV), we have f/f (C) > U2{0 for C e 
hiV,iV] . 

Consequently we have f/f (C) > ^^2(0 on R. 

Now, decrease 9 until one of the following situations happens. 

1. There exists a ^ > 0, such that C/f (C) = U2{£,)- In this case we have finished 
the proof. 

2. There exists a ^ > and Ci € such that one of the components of 
and U2 are equal there; and for all f G M, we have C/f (C) > C^2(C)- On applying 
the Maximum Principle on R for that component, we find J7f and U2 must be 
identical on that component. To fix ideas, we suppose that the component is the 
first component. Then Uf - U2 satisfies ([3125]), (p^ej) and ([3221) • Plugging wi = 
into p.26p we find that there is at least one Ce such that W2{£,g) = 0. Then by 
applying maximum principle to (I3.26p . we have W2(C) = for C € IR- Similarly we 
also find W3(C) = 0, C G R. We have then returned to case 1. 

Hence, in either situation, there exists a ^ > 0, such that 
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for all ^ e M. 

The nonexistence of the monotone traveling waves for comes from the fact 

that all its solutions are oscillatory for c < 2^/1 — ai. □ 

Concerning the wave solutions of system (11.21) we immediately have 

Corollary 14. Assume the conditions HI lll.6\) and H2a |i. 7p or conditions HI 
UTg]) . H2b fTTg]) and H3 UTgj) . T/ien for each c > 2^1 - ai and all t > system 
il.2\) has a unique monotonia traveling wave solution connecting the equilibrium 
(0, 1) with (1,0). 

The conclusion of the Corollary says the delay does not change the course of the 
traveling waves, but it may change the asymptotic behaviors of the wave solutions 
at +0O. 
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